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Abstract. We calculate the flux of neutrinos generated by the propagation of 
ultra-high energy nuclei over cosmological distances. The propagation takes into 
account the interactions with cosmic background radiations including the CMB 
and the most recent estimates of higher energy (infra-red, optical, and ultra violet) 
backgrounds. We assume that the composition of ultra-high energy cosmic rays 
(UHECRs) at the source is the same as the observed one at low energies. This 
assumption fits well the present data at the highest energies. We compare the 
cosmogenic neutrino flux from mixed composition sources to pure proton sources. 
We find that the neutrino flux in the mixed composition case has a high energy 
peak, mainly due to photopion production off CMB photons, of similar shape 
and amplitude to the proton case. At low energies both composition cases have 
significant neutrino flux with a peak around 10 14 ' 5 eV due to the higher energy 
backgrounds. The mixed composition case induces a higher flux of neutrinos 
at energies below 10 13 eV due to the neutron decay component that extends 
down to low energies. Detection of diffuse neutrino fluxes at ultra high energies 
can strongly constrain the source distribution of UHECR, whereas fluxes at lower 
energies could be used to constrain confinement of VHE and UHE cosmic rays if 
combined with composition analysis from cosmic ray experiments. 
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1. Introduction 

Mass composition measurements are key to solving the mystery of the origin of 
ultra- high energy cosmic rays (UHECRs). If UHECRs originate in extragalactic 
astrophysical accelerators, the observed composition of UHECRs should be primarily 
protons and light nuclei (see, e.g., P), while heavier nuclei could indicate a Galactic 
origin (see, e.g., [5]), and photon primaries could indicate top-down scenarios (see 
[3] for a review). Future composition measurements should be able to differentiate 
between alternative scenarios for the origin of UHECRs and to determine the transition 
between Galactic and extragalactic cosmic rays. 

Present studies of cosmic ray composition indicate a dominance of heavier nuclei 
(from 4 He to 56 Fe) around the knee region Ej (~ 10 16 -10 17 eV) followed by a tendency 
toward lighter nuclei around ~ 10 18 eV and above However, the composition 

of cosmic rays is notoriously difficult to determine at high energies. A number of 
airshower observables are composition dependent, such as the fraction of muons in 
the shower, the shower maximum, and the fluctuations of shower maximum. These 
observables require observatories with large statistics, control of systcmatics, and a 
large dynamic range in observed energies. In addition, comparisons of observations 
with model predictions depend on simulations of hadronic interactions in an energy 
and rapidity range that have not been tested by accelerator experiments. Current 
efforts in improving hadronic interaction models and the construction of the Pierre 
Auger Observatory [7] should lead to a much better understanding of UHECR 
composition. 

Another key observation that should become feasible in the next few years is 
the detection of cosmic neutrinos at high energies. UHECRs are of great interest as 
a source of high energy neutrinos that are almost guaranteed by their interactions 
with the cosmic background radiation. Cosmic rays at the highest energies produce 
pions off the cosmic microwave background (CMB) giving rise to the Greisen-Zatsepin- 
Kuzmin feature in the cosmic ray spectrum jHJ. The decay of the charged pions 
produced by the interactions with the CMB generates neutrinos IU] that are 
often called cosmogenic, GZK, or photopion neutrinos. Furthermore, it has been 
recently shown |12) . that the interactions of protons with the infra-red and optical 
intergalactic backgrounds produce a large amount of cosmogenic neutrinos at lower 
energies. 

Most predictions of the cosmogenic neutrino flux assume that the UHECR 
primaries are protons. More recently cosmogenic neutrino fluxes were calculated 
for other primaries such as pure 56 Fe ^3 ^] E3i 4 He, and 16 (TJ|. Here we 
calculate the cosmogenic neutrino flux for sources that inject a mixture of primaries 
with the same initial abundances as the observed Galactic cosmic rays and compare 
with the pure proton case. The mixed composition model was proposed in pQ and 
studied in detail in where a comparison between the predicted spectrum and 
composition for this model is contrasted with the pure proton assumption. We 
briefly review the mixed composition model and describe the adopted source evolution 
models in §2. In order to calculate the neutrino flux, we propagate the injected 
primaries and their daughter nuclei from cosmological distances to Earth calculating 
their interactions with cosmic backgrounds (§3). We simultaneously calculate the 
generated flux of neutrinos for composition at injection (pure proton or mixed) and 
different source evolution scenarios (§4). The generated flux of neutrinos for both 
composition hypotheses is normalized by comparing the predicted UHECR spectrum 
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with present observations. We discuss the implications of our results for ongoing and 
future experiments in §5 and conclude in §6. 

2. Mixed Composition UHECR and Source Evolution Models 

In order to compare the predicted cosmogenic neutrino flux for a mixed composition 
scenario with previous models, we study the neutrino production for two UHECR 
injection compositions: pure proton and mixed composition. The mixed nuclei case 
assumes that the injection composition at the source matches the abundances observed 
in Galactic cosmic rays at lower energies as in |17| . The source is assumed to emit a 
power law spectrum with spectral index a such that the number of nuclei i emitted 
in the energy range [E, E + dE] is: 

ni (E) = x i A^ 1 nE- a dE , (1) 

where A is the mass number of a given nucleus, k is a normalization constant, and Xi 
is the abundance of species i (given in |T7)). 

In the mixed composition model of yQ , the maximum energy at injection of each 
species is set to 

Emax.i — ZiE max ( H^j , (2) 

i.e., the maximun energy is proportional to the charge number Z of a given nucleus. 
This assumption is reasonable if the maximun energy at the source is limited by 
the confinment of particles. The evolution of the relative abundances under these 
assumptions as a function of the injection spectral index is displayed in Fig.^ For the 
neutrino flux calculation, we first set the maximum proton energy E max ( l H) = 10 20 ' 5 
eV for the pure proton case and use the same input E max ( l H) in eq. El to set the 
mixed composition case, which gives a maximum energy for iron of ~ 10 22 eV. We 
assume an exponential cut-off at the source above E max (Z). We discuss the effect of 
changes in E max to the expected neutrino flux below. 




beta (logarithmic index) 



Figure 1. Evolution of the integrated relative abundances at the source 
normalized to protons as a function of the injection spectral index between 

Emax/ 10 and Ernax- 



The predicted neutrino flux is strongly dependent on the choice of source redshift 
evolution. We consider four source evolution models in our predictions. First, we 
assume no evolution with redshift (hereafter called uniform distribution). Second, we 
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consider a redshift dependence proportional to the old estimate of the star formation 
rate (SFR) of [TH] (hereafter oSFR), which evolves as (l + z) n for z < 1.9, followed by a 
constant (2.9)™ between 1.9 < z < 2.7, and an exponential cutoff (2.9) n exp(l — z/2.7) 
for z > 2.7 and with n — 3 as assumed by most previous cosmogenic neutrino flux 
calculations. Our third assumption is based on a more recent estimate of the SFR 
evolution (hereafter new SFR or nSFR) that can be deduced from pD] and evolves as 
(1 + z) 3 for z < 1.3, followed by a constant rate between 1.3 < z < 6, followed by 
a sharp cut-off. Finally, we use a stronger source evolution (hereafter called strong 
evolution) favored by the recent infra-red survey of the Spitzer telescope |2] ■ We use 
the following parametrization of their model: (1 + z) 4 for z < 1, followed by a constant 
rate between 1 < z < 6, followed by a sharp cut-off. The two latter source evolution 
hypothesis were recently used in |12j . 

3. Interactions of protons and nuclei with photon backgrounds and 
neutrino emission 

In the following sections we consider the interaction of protons and nuclei with the 
CMB and the infra-red, optical and ultraviolet backgrounds (hereafter we group these 
three backgrounds under IR/Opt/UV for short). The effect of IR/Opt/UV photons 
on the propagated UHECR spectrum of pure proton sources is negligible, however, as 
shown in ^21> these additional backgrounds have a significant effect on the neutrino 
flux associated with UHECRs. To model the IR/Opt/UV backgrounds and their 
cosmological evolution, we use the latest estimate of |22I which is based the earlier 
work in Refs. [51] and updated with recent data on history of the star formation 
rate and the evolution of galaxy luminosity functions. We use IR/Opt/UV calculated 
at 26 different redshifts (Az = 0.2) between and 5. The differential density of the 
IR/Opt/UV is shown on Fig. [2^. In general, the IR/Opt/UV has a much milder 
cosmological evolution when compared to the CMB. However, the evolution of the 
photon background in the UV range is notable for redshifts between 1 and 0. The UV 
background gets fainter at lower redshifts due to the aging of the stellar population 
which results in the death of the stars with the shortest lifetimes which produce all 
of the UV emission [21]. The basic reason is that most galaxies were forming stars at 
much higher rates at z=l than they are today (e.g., [231) ■ The UV background has a 
sharp cut-off above 13.6 eV due to the Lyman limit absorption from HI observed in 
galaxies PS]. In the following, we assume this edge to be black, i.e., we assume that 
there are no photons above 13.6 eV. 

Protons and nuclei propagating in the extragalactic medium interact with CMB 
and IR/Opt/UV background photons. These interactions produce features in the 
propagated UHECR spectrum such as the GZK cutoff [BJ and their decay products 
generate the cosmogenic neutrino flux. The proton mean free paths for photo-pion 
interactions with the CMB and IR/Opt/UV background at different redshifts are 
displayed on Fig. [2p. The evolution of the mean free path with redshift is strong in 
the case of the CMB, which implies that the dominant background strongly depends 
on the source evolution assumed. It is also important to note that the photo-pion 
production off IR/Opt/UV photons competes with the pair production process off 
CMB photons which is also evolving with redshift. Due to the slow evolution of the 
density of the IR/Opt/UV background with the energy, the mean free path above the 
interaction threshold (determined by the UV break of the background) evolves slowly. 
A sharper decrease of the mean free path is visible at higher energies when protons 
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Figure 2. Left: Evolution of the differential density of the IR/Opt/UV 
background for redshifts between and 5 1221 , Right: Evolution of the proton 
mean free path for photo-pion production with the IR/Opt/UV and the CMB at 
redshifts between and 5 

start to interaction with photons of the far IR bump. At very high energies (above 
~ 5 x 10 19 (1 + z) eV) the CMB contribution starts to dominate the mean free path 
evolution and the effect of the IR/Opt/UV photons becomes negligible. 

Fig. shows the contribution of the different processes with the different 
backgrounds to the total attenuation length of protons at z = 0. This plot clearly 
shows how the photopion production process off the IR/Opt/UV background has little 
effect on the predicted UHE proton spectra. Indeed, as soon as the pair production 
threshold with the CMB is reached, the photopion production with IR/Opt/UV 
photons is completely subdominant and can be neglected in the calculation of UHECR 
with pure proton sources. 

The interactions experienced by nuclei with photon backgrounds are different from 
the proton case. Pair production (for which we use the mass and charge scaling given in 
\'27\) results in a decrease of the Lorentz factor of the UHE nucleus, whereas at higher 
energies, photodisintegration (also called photoerosion) processes lead to the ejection 
of one or several nuclcons from the nucleus. Different photoerosion processes become 
dominant in the total interaction cross section at different energies [2E1- The lowest 
energy disintegration process is the Giant Dipole Resonance (GDR) which results in 
the emission of one or two nucleons and a particles. The GDR process is the most 
relevant as it has the highest cross section with thresholds between 10 and 20 MeV 
for all nuclei. For nuclei with mass A > 9, we use the theoretically calculated GDR 
cross sections presented in |29) . which take into account all the individual reaction 
channels (n, p, 2n, 2p, np, a,...) and are in somewhat better agreement with data 
than previous treatments. For nuclei with A < 9, we use the phenomcnological fits 
to the data provided by [27j. Around 30 MeV in the nucleus rest frame and up 
to the photopion production threshold, the quasi-deuteron (QD) process becomes 
comparable to the GDR and dominates the total cross section at higher energies. The 
photopion production (or baryonic resonances (BR)) of nuclei becomes relevant above 
150 MeV in the nuclei rest frame (e.g., ~ 5 x 10 21 eV in the lab frame for iron nuclei 
interacting with the CMB), and we use the parametrisation given in |27j where the 
cross section in this energy range is proportional to the mass of the nucleus (nuclear 
shadowing effects are expected to break this scaling above 1 GeV). The reference for 
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this scaling is the deuteron photoabsorption cross section which is known in great 
detail. It is important to note that photopion cross sections for nuclei are different 
from the free nucleon case. In particular, in nuclei the baryonic resonances heavier 
than the first A resonance are far less pronounced than for nucleons and the cross 
sections are not simply derived from the free nucleon case. We also follow |57| for the 
treatment of the nucleon multiplicities, the energy losses and absorption probability of 
the produced pion in the parent nucleus (see below). Finally, above 1 GeV, the total 
cross section is dominated by the photofragmentation process which fragments nuclei 
into individual nucleons or low mass nuclei. For the choices of E max in the present 
work, the photofragmentation process is negligible. 




log(y) log(y) 

Figure 3. Loft: Evolution of the attenuation length of protons as a function of 
the energy at z = 0. The contribution of pair production and pion production off 
the CMB and IR/Opt/UV are separated. Right: Evolution of the iron nucleus 
mean free path for the different photoerosion processes and interactions with the 
CMB and IR/Opt/UV photons at z = 

The contribution of the different photoerosion processes and the different 
backgrounds to the total mean free path for iron nuclei are displayed in Fig. |3t>. The 
photoerosion is dominated by the GDR process through most of the Lorentz factor 
range. The baryonic resonances begin to dominate only above 10 215 eV where the 
effect of the GDR starts to decrease. The BR process off the IR/Opt/UV background 
does not affect the propagated cosmic ray spectrum, however, as we discuss below, its 
contribution to neutrino production is important. 

Fig. 2t shows the contribution of pair production and photoerosion processes to 
the total attenuation length of iron nuclei. Photoerosion processes dominate through 
most of the energy range and the effect of pair production is small at low redshifts. 
Although the competition between pair production off the CMB and photoerosion 
processes with IR/Opt/UV photons depends on the redshift (e.g., at high redshifts 
pair production increases due to the strong evolution of the CMB), the propagation 
of nuclei is mainly dominated by photoerosion processes. A comparison between the 
attenuation lengths of different species is displayed in Fig.^J). The figure shows what is 
known since |28j , that the attenuation length of low mass nuclei are smaller than that 
of protons and heavy nuclei and, as a consequence, light nuclei should not contribute 
as significantly at the high energy end of the spectrum. Furthermore, iron nuclei 
have larger or similar attenuation lengths to protons up to 3 x 10 20 eV. However, the 
energy loss processes are different for protons and nuclei and the sole comparison of 
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attenuation lengths can be misleading. Most of the energy losses of nuclei result in 
nucleon ejection, thus, unlike protons, a given nucleus does not remain on "the same 
attenuation length curve" during its propagation. Therefore, in a mixed composition 
model, the iron component starts to drop around 3 — 4 x 10 19 eV due to the reduction 
of the GDR mean free path off far IR photons, whereas the proton component starts 
to increase. 
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Figure 4. Left: Evolution of the attenuation length of iron as a function of the 
energy at z = 0. The contribution of pair production and photoerosion processes 
off the CMB and IR/Opt/UV photons are separated. Right: Comparison of the 
attenuation length of different nuclei at z = 0. 

We use a Monte-Carlo code to propagate nuclei from the source to Earth as 
described in detail in Q]. The neutrinos produced during the propagation of each 
particle are recorded with their flavor, energy, and production redshift. Neutrinos are 
produced by the decay of pions and of secondary neutrons. In the case of secondary 
neutron decay, n — > p + e~ + 77 e , the energies of the outgoing particles are calculated 
with a three body decay algorithm. 

Neutrinos are also produced by the photopion production of protons and neutrons. 
This process has been treated in great detail in jHO], where all the relevant baryonic 
resonances and possible mesons and multi-pions channels were taken into account using 
|31| . In this work, we use a simpler treatment assuming that the total cross section is 
dominated by the emission of single pions and use a classical model of the kinematics 
of the delta resonance process: N + j-^A^N' + tt (where N is a nucleon). For the 
energy range we consider, the delta resonance dominates multi-pion production and 
most neutrinos are produced close to threshold [32], ^UJ- As we discuss below, the 
shape of the neutrino spectrum we obtain closely follows the more detailed treatment 
of PHI (see, e.g., Fig. HJ). Photopion production through the delta resonance has a 
1/3 probability of isospin flip of the incoming nucleon, and each isospin flip leads to 
the production of three neutrinos. For example, in the case of proton interactions 
producing tt + , p + 7 — > ir + + n, the ir + decay generates one Vp,, one v e + and one 
F M . The charged pion decay, ir + — > ji^ -\- , is calculated using the two body decay 
algorithm, while for the muon decay the energies are calculated using the three body 
decay algorithm: /i + — > e + + v e + 

We use the model above to calculate the neutrino production from interactions of 
primary and secondary nucleons and nuclei with the CMB and IR/Opt/UV photons. 
In the case of nuclei propagation, neutrinos can also be produced via the photopion 
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production of nuclei. The treatment of this component is more uncertain complicated 
by pion interactions within the nucleus. To a good approximation, the pion production 
can be treated as the production of free nucleons, but the energy losses of the 
pion inside the nucleus have to be taken into account. The transfer of the initial 
kinetic energy of the pion to the spectator nucleons is partly responsible for the high 
multiplicity of ejected nucleons from this process |22]. Furthermore, the produced 
pion can also be absorbed by a pair of nucleons before leaving the nucleus. To 
model these different effects, we follow 27i and use an absorption probability for 
the produced pion which is proportional to the nucleus radius P a bs — 0.13A 1 ' 3 . 
This approach reproduces well the available data for a small number of nuclei. If 
the pion is not absorbed, we remove 40 MeV of kinetic energy in the nucleus rest 
frame per participant nucleon before calculating the energy of the produced neutrinos. 
Although this simple treatment allows a reasonable estimate of the neutrino flux due 
to photopion production, a more sophisticated treatment would be useful for more 
precise calculations in the future. 

4. Neutrino fluxes for cosmological distributions of sources 

Fit to the observed UHECR spectrum 

We used our monte-carlo simulation to calculate the differential flux of neutrinos for 
different cosmological distribution of sources. We propagate samples of 10 8 particles 
between E min = 10 16 eV and E max = Z x 10 20 5 eV. The particles for sources 
between z — 10~ 5 and z = 6 are followed down to 10 15 eV where the interaction 
probability becomes negligible. We inject protons and a mixed composition of nuclei 
and propagate the particles assuming different source evolutions as described above. 
For each composition hypothesis, the nominal value of the injection spectral index, 
0, is chosen by requiring a good fit between the simulated and the observed UHECR 
spectrum. Since the neutrino production is calculated within the same Monte-Carlo 
code, the normalization of the simulated spectra to the observed one provides the 
normalization for the neutrino fluxes. We normalize the UHECR flux to the HiRes 
and Auger [33] spectra, which are a factor of 1.8 below the AG AS A 35 : spectrum 
at 10 19 eV. 

The propagated UHECR spectra that best fit the data for both composition 
hypotheses and different source evolutions are displayed on Fig. As shown in 
previous studies (e.g., 123 EZ|)> the generated proton spectra and the presence of 
a pair production dip are only mildly dependent on the source evolution hypothesis. 
The amplitude and the energy of the pair production dip minimum (around 10 18 ' 7 eV) 
only slightly depend on the evolution. However, the beginning of the dip is determined 
by the transition between the adiabatic and the pair production losses, which is more 
sensitive to the density of sources at high redshifts [SSj- We assume that the galactic 
component compensates for the difference between the predicted flux and the observed 
one at energies below the dip. We discuss the implications of the different galactic 
to extragalactic transition and source evolution models for the neutrino flux below 
(§4-2). 

In the mixed composition case the transition from galactic to extragalactic 
components ends at the ankle. In Fig. |3Jd we show the results for three source 
evolution models (uniform, oSFR, strong). The transition point and the spectrum 
above 10 1 eV is quite insensitive to the source distribution. At energies below the 
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ankle, the galactic fraction that completes the total observed spectrum depends on 
the source evolution, as in the case of pure proton models below 10 18 eV. Although 
the shape of the predicted spectra and their implications for composition are not 
dependent on the source evolution, the spectral indices required to fit the data are 
harder for stronger source evolution. For the pure proton model (3 moves from 2.6 in 
the uniform evolution case to 2.4 for the strong evolution case (/?=2.5 for oSFR and 
nSFR) [2E1 CHU OH EH \U\ , while in the mixed composition scenario (3 varies from 2.3 
for the uniform evolution case to 2.1 for the strong evolution case ((3=2.2 for oSFR 
and nSFR). 
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Figure 5. Left: Propagated spectra obtained assuming a pure proton 
composition and four different source evolution hypotheses compared to the HiRcs 
monocular data at ultra high energy. Right: Same as left but assuming a mixed 
composition and three source evolution hypotheses. 



4-2. Neutrino fluxes for pure proton sources 

For each composition and source evolution hypothesis our propagation code calculates 
a normalized flux of neutrinos produced during the propagation. We assume that 
the produced neutrinos only lose energy due the expansion of the Universe and 
we neglect their interactions with other particles. Before estimating the neutrino 
flux for a given UHECR spectrum, we checked that our simplified treatment of the 
neutrino production, which uses only the single pion production, gives an accurate 
estimate of the neutrino flux. We calculated the neutrino flux with our code under 
the same assumptions of [3D], i.e., only interactions with CMB photons, the oSFR 
source evolution model, a injection spectral index [3 = 2.0, and a maximum energy at 
the source of E max = 10 215 eV. A comparison between the two studies is displayed 
in Fig. El where the flux of z/ M 's produced by photopion interactions is plotted. 
The differences between the two calculation are very small never exceeding ~ 10%. 
This comparison shows that our approximation is accurate enough for a comparison 
between neutrinos from different injected compositions, especially considering the 
other uncertainties in the assumptions (e.g., spectral index, E max , source evolution, 
and magnetic fields). 

Neutrino fluxes for a pure proton composition and four different source evolutions 
are displayed on Fig. [7Ji where it is clear that the source evolution is a critical 
parameter. At high energies, ~ 10 18 eV, the neutrino flux for a uniform distribution 
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Figure 6. Comparison of the muon neutrino flux calculated in the present 
work and in 1201 for the oSFR source evolution, a spectral index f3 = 2.0, and 
E max = 10 21 - 5 eV. 



of sources is almost one order of magnitude below the other hypotheses 021 QUI 
If the sources of UHECRs are distributed uniformly in redshift, the observation of 
cosmogenic neutrinos will be very challenging unless the sources can reach extremely 
high maximum energies |48j . For the three other source evolutions, the neutrino fluxes 
are comparable: the oSFR source evolution gives a stronger weight to high redshift 
sources, but in the strong evolution case the harder spectral index counterbalances 
the redshift effect. 




Figure 7. Left: Total neutrino flux from pure proton sources and four source 
evolution hypotheses. Right: Contribution of the different backgrounds to the 
neutrino flux for a strong source evolution hypothesis. The photopion interactions 
with CMB (dashed line) and IR/Opt/UV (dotted line) are shown and the 
difference between the sum of the two contributions and the total is the neutron 
decay component. 

The contribution of the CMB and IR/Opt/UV backgrounds is detailed for the 
strong evolution case in Fig. [TJz>. At low energies, the neutrino flux is dominated by 
the contribution of the IR/Opt/UV backgrounds ^2] an( i the flux is much higher 
than at higher energies. Although the interaction probability with the IR/Opt/UV 
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backgrounds is much lower than with the CMB, the number of particles that are able 
to interact is much higher due to the steep injection spectra provided by acceleration 
mechanisms in astrophysical sources. The CMB contribution generates a peak at 
~ 10 17 - 6 eV, while the peak at ~ 10 14 5 eV is due to the IR/Opt/UV contribution. The 
position of the peaks depends on the combination of the evolution of the interaction 
probability and the injection spectrum. In Fig. above the threshold, the mean 
path of protons (thus, the interaction probability) evolves very slowly with energy. 
As the injection spectra are steep (i.e., more particles are injected at low energies), 
the number of protons interacting by photopion production with IR/Opt/UV photons 
is continuously decreasing above the threshold (~ 10 16 eV at z = 1), the neutrino 
peak is then due to interactions close to the threshold corresponding to neutrino 
energies around 10 14 ' 5 eV. The interaction probability dramatically increases at higher 
energies, above the interaction threshold with CMB photons, resulting in a break of 
the evolution of the flux and a high energy neutrino peak. The peak at high energy 
is mainly due to neutrinos produced close to the threshold, therefore, the energy of 
the peak is only mildly dependent on the maximum energy at the sources. Due to 
the slow evolution of the IR/Opt/UV background with redshift, the effect of source 
evolution on the neutrino flux is milder than for the CMB neutrinos, but the influence 
of the spectral index is higher. The neutrino flux at lower energies are higher for 
steeper injection spectra, which explains why the oSFR and nSFR evolution scenarios 
generate higher flux than the strong evolution shown in 01. 




12 13 14 15 16 17 16 M 20 14 15 16 17 16 19 20 
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Figure 8. Left: Contribution of different neutrino flavors to the total flux for 
strong source evolution. Plotted are P e 's from neutron decay and (+ u^) and 
Ve (+ v e ) from pion production. Right: Contributions from different redshifts to 
the neutrino flux for strong source evolution. 



The generated flux of different flavors of neutrinos is displayed on Fig. [SJi for a 
strong source evolution hypothesis. (Note that the observed fluxes will be affected by 
mixing.) Each component (ia, + and v e + v e from pion production and v e from 
neutron /3-decay) exhibits a double peak structure due to the contribution from the 
CMB at high energies and the IR/Opt/UV backgrounds at low energies. The neutron 
decay component peaks are shifted to lower energies due to the kinematics of the 
process. The contribution at low energies coming from neutron decay after proton 
interactions with UV photons induces the third peak in the total neutrino flux around 
10 12 - 5 eV. 
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The contribution from different redshifts to the production of cosmogenic 
neutrinos is displayed on Fig. [SJd for the strong source evolution hypothesis. It is clear 
that the contributions from different redshifts depend on the chosen background since 
the backgrounds evolve differently with redshift. At high energies, the contribution 
from high redshifts (z > 1.5) is important due to the strong evolution of the CMB and 
the sources. The peaks generated at different redshifts move in energy mainly due to 
the evolution of the temperature of the background that allow lower energy protons 
to interact at higher redshifts. The strong evolution of the neutrino production with 
redshift explains the slow increase of the neutrino flux below the high energy peak. In 
the high energy peak region (around 10 17 ' 6 eV) most of the neutrinos are produced 
at reshifts above 1.5. In the case of the IR/Opt/UV component, the dominant 
contribution at the low energy peak comes from intermediate redshifts (between 0.5 
and 1.5) as the slow evolution of the backgrounds cannot compensate for the decrease 
of the number of sources at high redshifts. For the three redshift ranges above z = 0.5, 
the neutrino flux exhibits peaks at similar energies (the differences are due to adiabatic 
loses rather than the evolution of the background which is mild in the UV range for 
z >1). The low redshift contribution (z < 0.5) is quite different, as the flux peak from 
starlight is shifted to higher energies, due to the low-redshift drop in UV emission 
from young stars evolving off the main sequence. 



4-3. Neutrino fluxes for a mixed composition 

In the case of a mixed composition, the expected neutrino flux is shown in Fig. for 
the different source evolution hypotheses. At high energies, the flux is very similar 
to the pure proton case, which is not surprising since the mixed composition models 
with j3 between 2.1 and 2.3 are proton dominated (see Fig. The flux is much 
higher in the strong evolution case compared to a uniform distribution. There is very 
little difference between the strong and oSTR evolution cases (for simplicity, we omit 
the nSFR). The intermediate energy peak appears around 10 14 ' 5 eV and the neutron 
decay peak is just below 10 12 eV. 




Figure 9. Left: Total neutrino fluxes corresponding a mixed composition and 
three source evolution hypotheses (the nSFR, hypothesis is omitted). Right: 
Contribution of the different species to the neutrino flux for a strong source 
evolution hypothesis 



In the case of a strong evolution, the contribution of protons, He, and Fe nuclei 
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to the total flux are displayed on Fig. 03. Above 10 15 eV, the main contribution is due 
to protons (~ 65%), the second contribution comes from He, and Fe nuclei contribute 
only of a few percent over the whole energy range. Due to the harder spectral index, 
the intermediate energies peak for the proton component in Fig. [SJd is slightly shifted 
towards higher energies (to ~ 10 15 eV), when compared to the corresponding peak in 
the pure proton case. The intermediate energy peak for the total flux from all the 
species is again around ~ 10 14 ' 5 eV due to the contribution of the other species which 
become dominant for energies below the proton peak. The neutron decay peak of the 
proton component is invisible on the total flux as it is completely dominated by the 
nuclei contribution in this energy range. It is important to note however that the 
contribution of the different species is strongly dependent on the relative abundances 
assumed at the source. In our model, protons contribute ~ 50%, helium ~ 30%, and 
Fe ~ 5% at the source for a spectral index (3 = 2.1. 

To show the relative contributions to the neutrino flux from different species, we 
show in Fig. HOfa the generated fluxes in arbitrary units assuming the same spectral 
index and similar integrated abundances at the sources between 10 16 and Z10 20 5 eV 
for all the species. It is clear that, except for protons, there is very little difference 
between the contributions from different species under the assumption of a maximum 
energy scaling with Z. As mentioned in JH], the acceleration of nuclei above 10 20 eV 
is not required to fit the UHECR spectrum if protons are accelerated above ~ 10 20 
eV, therefore, the contribution of nuclei at high energies is uncertain and strongly 
dependent on the maximum energy reached at the source. 

In Fig. 110b we show the contribution of the different production processes for 
He nuclei which is qualitatively similar to those of heavier nuclei. The high energy 
neutrino production for nuclei is dominated by the pion production of secondary 
nucleons off CMB photons, which requires nuclei to be accelerated above ~ A x 
5 x 10 19 eV/ (1 + z) corresponding to E max ~ 10 21 eV for Fe nuclei. For the maximum 
energy we assumed (see §2), the flux from direct photopion production of nuclei is low 
due to the fact that this process only becomes dominant above ~ A x 1.5 X 10 20 eV 
and that the produced pion can be absorbed before producing neutrinos (the relative 
contribution from this process for iron is even lower than in the case of helium). If 
nuclei are not accelerated to the highest energies the flux is lower than our estimate, 
while the flux could be slightly higher if the proton abundance is higher than in our 
model. 

In the intermediate peak region around 10 14 ' 5 eV, the neutrino flux from nuclei 
originate mainly from direct pion production and neutron decay from secondary nuclei. 
Direct photopion production is the dominant process in the peak region, which may 
appear surprising since the GDR interaction probability (which is responsible of the 
nucleon emission) is higher than the photopion probability. However, for the photopion 
production of secondary nucleons, the neutrino emission requires first the emission of 
a nucleon and a subsequent photopion interaction of the emitted nucleon. At low 
energies, the realization of both of these requirements becomes less probable than a 
direct photopion interaction (despite the pion absorption probability) , therefore, direct 
photopion dominates. The neutrino flux produced by secondary nucleon photopion 
production drops faster at low energies than in the free nucleon case (see the proton 
component on Fig. HUb). 

At the lower energies, below 10 14 eV, the neutrino flux from nuclei is dominated 
by the contribution of secondary neutron /3-decays. This component is far more 
important than in the case of primary protons as the interaction probability of nuclei 
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Figure 10. Left: Comparison of the neutrino production of different species 
assuming the same integrated abundances for all species between 10 16 eV and 
Z 10 20 ' 5 eV, a strong source evolution scenario, and /3 = 2.1. Right: Contribution 
of the different processes of neutrino production for the He generated neutrino 
flux in the left panel. 



(via the GDR process) is much higher. Unlike the secondary nucleon photopion 
component, the flux keeps increasing at lower energies as a subsequent interaction 
of the ejected nucleon is not necessary. The CMB contribution peaks around 10 14 
eV but is overwhelmed by the IR/Opt/UV background contribution in the case of 
He nuclei. The shift of the neutron decay peak to lower energies compared with the 
pure proton case can be explained by the fact that the energy per nucleon threshold 
of the GDR process is lower than for the photopion process. Neutrons with energies 
below 10 16 eV can be ejected in contrast with the pure proton case where neutrons are 
produced following a photopion interaction of a proton, which does not occurs below 
10 16 eV even at high redshifts. 

Finally, when comparing the contribution of the different species it is important to 
note the sensitivity of the generated flux to the spectral index assumed. Indeed, as the 
threshold energies for the different processes scale approximately with the mass of the 
nuclei (i.e., they occur at approximately the same Lorentz factor for all species), lower 
mass nuclei have lower energy thresholds. On the other hand, the cross section and the 
number of nucleons that can be ejected are higher for higher mass nuclei. The steep 
spectral indices in astrophysical accelerators together with lower energy thresholds 
favor the neutrino production of low mass nuclei, which is why the neutrino production 
is higher for He nuclei. The evolution of the cross sections and of the number of 
nucleons with the mass tend to counterbalance partially this effect especially if the 
spectral index is hard. For (3 — 2.1 the difference between the neutrino production 
from different nuclei is relatively small, while for softer spectral indices the lower mass 
nuclei become even more dominant over the whole energy range. 

4-4- Comparison between the two composition hypotheses 

The comparison between the pure proton and the mixed composition is displayed 
on Fig. ^2 for a strong source evolution hypothesis. For the other choices of source 
evolution, the comparison is qualitatively similar. At high energies, where detection 
is feasible, the fluxes are comparable. Changes in the detailed composition at the 
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source or in the maximum energy reachable for the nuclei can slightly increase or 
decrease the predictions, but for a proton dominated mixed composition the flux at 
the highest energies is very close to the pure proton model. The similarity between 
the fluxes arises from the fact that the lower contribution of nuclei to the neutrino 
flux is compensated by the harder spectral index required to fit the UHECR data. 
Conversely, at lower energies, the higher fluxes expected for the pure proton scenario 
are mainly due to the softer spectral index that give a higher luminosity at low 
energies. In the intermediate energy range, detectability is limited by the atmosphere 
neutrino background for energies below 10 15 eV. For energies between 10 15 eV and 
10 17 eV, the pure proton case increases the chances of future detectability. However, 
the comparison between the two compositions is sensitive to a few assumptions as we 
discuss below. 




log 10 E(eV) 

Figure 11. Predicted cosmogenic neutrino flux for pure proton and mixed 
composition in the case of the strong source evolution compared with the 
sensitivities of Auger (u T ) and ANITA at high energy and the limit of Ice Cube 
(for the Vp detection channel only) after three years of observation with 90 C.L. 
(assuming neutrino spectra tx E~ 2 and E ) at low energy, estimates of the 
atmospheric neutrino flux and atmospheric neutrinos due to charmed interactions 
are also displayed. 



5. Discussion 

In the previous sections, we calculated neutrino fluxes for mixed and pure proton 
compositions for different source evolution hypotheses. In this work, we made some 
assumptions prior to the calculation that we further discuss below as they raise 
some uncertainties on the expected cosmogenic neutrino fluxes. We also discuss 
the IR/Opt/UV background and the possible detection of the predicted cosmogenic 
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neutrino fluxes. 

5.1. Change of the proton injection spectral index at low energy 

In the pure proton case, we assumed that the injection spectrum was constant down 
to low energies. This assumption affects the transition between the galactic and 
extragalactic cosmic ray components. In particular, the extragalactic proton flux 
obtained around 10 17 eV is too high to match the low proton abundance observed 
by KASCADE 4 . A change of the injection spectral index around 10 18 , making it 
harder at low energies, can alleviate the problem with the Kascade data and prevent 
the sources from being too luminous as proposed by |39) . In this model, the transition 
between heavy galactic and proton extragalactic components is extremely steep as 
discussed in |16) . The consequence of such a change in the spectral index for the 
neutrino flux in the strong evolution case is to bring the pure proton case down to 
mixed composition case, i.e., a flux well below the one displayed in Fig.0 

5. 2. Neutrino fluxes and extragalactic magnetic fields 

Another way to solve the problem of the proton abundance at 10 17 eV is to invoke 
a magnetic horizon effect 0H] that prevents lower energy extragalactic protons 
from reaching the Earth. The effect on the steepness of the galactic to extragalactic 
transition is similar to the proposal discussed above (§5.1), but without the need to 
change the spectral index. The neutrino fluxes we calculated would not be significantly 
changed as the magnetic fields would only have an effect on the charged particles. 

Strong extragalactic magnetic fields can change the propagated UHECR spectrum 
and may challenge the fits in §4, specially for mixed composition models The 
fairly weak observational and theoretical constraints on the intensity and configuration 
of extragalactic magnetic fields make a precise estimate of their possible effects very 
difficult (see, e.g., The effect of magnetic fields should be most significant at 

lower energies, where we assumed that the particles could escape the source for both 
composition hypotheses down to 10 16 eV. This assumption is somewhat extreme even 
for pure protons. Astrophysical sources are expected to be surrounded by magnetic 
fields and radiation fields that can be quite strong, for example in galaxy clusters. 
Charged particles can be confined and the diffuse neutrino flux in the low energy 
range could be significantly higher than the fluxes we calculated above. 

5.3. Infrared, Optical, and Ultraviolet Backgrounds 

In order to calculate intergalactic IR photon fluxes and densities and their evolution 
over time (or redshift), the authors of Ref. |22 adopted the approach of using an 
empirically based method of calculating the infrared background radiation based 
on (1) the luminosity dependent galaxy spectral energy distributions obtained from 
observations, (2) observationally derived galaxy luminosity distribution functions and 
(3) the latest redshift-dependent luminosity evolution functions, primarily based on 
recent data from the Spitzer infrared space telescope. As is shown in Figure 6 of that 
paper, the predicted infrared flux distribution as a function of energy hugs the data 
from the Cosmic Background Explorer within 1 a, the lower limits from galaxy counts, 
and the upper limit from TeV 7-ray constraints. The estimated photon backgrounds 
emitted by stars and warmed dust grains in galaxies are now well measured from the 
current epoch back to a redshift of 1 . The uncertainly in the far infrared energy range 
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is ^20% and that in the mid-IR range is ~30%. The two different models in Ref. 22 
differ in the optical-UV range by ~ 50%. These are a fair measures of the observational 
and calculational uncertainties in these energy ranges. However, at redshifts greater 
than 1 the uncertainty in the integrated galactic photon emission could be significantly 
larger, as is the uncertainty in the ratio of far-infrared to UV emission from dust 
absorption and reradiation at these redshifts. However, the strongest feature in the 
background photon spectrum, viz., the sharp cutoff below the Lyman limit responsible 
for the intermediate peak, has a more solid theoretical and observational foundation. 
None of the dozen or more galaxies imaged above this limit of 13.6 eV energy has any 
detected continuum radiation, even at the level of a few percent of the emission at 
energies below this cutoff (EH). Although, the uncertainties on the backgrounds are 
higher at high redshifts, the evolution of the IR-UV backgrounds with redshift remains 
milder than in the case of the CMB. Therefore the contribution of high redshifts to 
the neutrino fluxes would anyway not be dominant even with a significantly higher 
background. As the IR-UV backgrounds are well constrained at low and intermediate 
redshifts [z < 1.5), we can conclude that the knowledge these backgrounds is good 
enough not to be a major source of uncertainties for our calculations. 

5.4- Detection of the cosmogenic neutrino fluxes 

The neutrino fluxes that we obtained assuming a strong source evolution for pure 
proton and proton dominated mixed composition models are displayed on Fig. 1111 In 
this figure we use the AGASA normalization which gives neutrino fluxes about 80% 
higher than our earlier estimates. This 80% difference between HiRes and AGASA 
normalization gives a fair estimate of the uncertainties triggered by the choice of the 
normalization. In the figure we also show the sensitivity of Ice Cube jI!5J for three years 
of observations and two different assumptions for the spectral index of the neutrino 
flux: = 2.0 and /3=1.5, the sensitivity of ANITA |F0] for a 45 days flight time, 
and the sensitivity of Auger to z^ T 's ^ETXj assuming maximum mixing and a sensitivity 
ranging from one event a year to one a decade. 

For the strong source evolution hypothesis, the flux that we calculated is close to 
the ANITA sensitivity and would be detectable if the maximum energy of acceleration 
is high enough, as it can be seen in the case of the mixed composition for E max = 
10 21 ' 5 eV. The typical case treated in most neutrino studies assumes a pure proton 
composition with a spectral index (3 = 2.0, the oSFR source evolution, and a high 
value of the maximum energy (typically 10 21 ' 5 eV), which generally give higher fluxes 
than the ones we calculated here and slightly above the ANITA sensitivity. However, 
this typical scenario with a hard injection spectral index does not provide the best 
fit to the UHECR spectrum, especially below 10 19 eV. It thus requires a Galactic 
component extending up to high energies, with typically equal flux of the galactic and 
extragalactic components at the ankle. The differences in the neutrino fluxes obtained 
in each case appear too weak to allow one constraining the nature of the ankle with 
the sole observation of high energy neutrinos. 

Since the predicted fluxes in the case of a mild source evolution are unfortunately 
too low to be observed by present detectors for reasonable assumptions on the 
maximum energies, the observation of any UHE neutrinos by Auger or ANITA would 
put very severe constrains on the source evolution |42j . which would then have to be 
strong or very strong (and coupled to a high value of the maximum energy). Note 
that the distribution of plausible sites of UHECR sources such as active galactic nuclei, 
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gamma-ray bursts, and galaxy clusters, evolve strongly with redshift, therefore one 
can conclude that, despite the high value of maximum energy required, the possibility 
of detecting neutrinos at high energy is not unlikely. 

In the low energy range, the situation is more complicated since the amplitude 
of the cosmogenic diffuse flux depends on the confinement efficiency at the source. 
However, the fluxes predicted assuming that the particles propagate freely in the 
intergalactic space are reasonably high. In Fig. 1111 we compare the predicted fluxes 
with the 90 % confidence level limits that Ice Cube should provide after three years 
of observations, for two different assumptions on the spectral index of the neutrino 
flux, (3 — 2.0 and /3=1.5, which correspond roughly to the proton and the mixed 
cases, respectively. Note that the limits shown only consider the channel, whereas 
Ice Cube should also be sensitive to electromagnetic showers from v r and v e |52| . 
Although the fluxes at lower energies are not as close to the instruments sensitivity 
as they are at high energies, and even more so for a mixed composition, they could 
be greatly enhanced if particles would be confined for a sufficiently long time around 
the sources. This is actually easier to achieve in the case of nuclei. Therefore, if Ice 
Cube observes a diffuse flux, it may be due to the confinement of nuclei in VHE and 
UHECR sources. 

Due to its large sensitivity, Ice Cube may be able to constrain the acceleration 
and the confinement of cosmic rays in sources even if cosmic ray nuclei are present. 
At low energies, the detection of a diffuse flux is complicated by the presence 
of the atmospheric neutrino flux and the uncertain flux from charm decay (see 
Fig. 1 1 1 ft . Methods to deconvolve atmospheric neutrinos from astrophysical diffuse 
flux were proposed in and may be effective in pulling a diffuse flux in different 
detection channels. Neutrino point source detections would be even more illuminating 
allowing for a significant improvement in studying UHECR confinement and sources 
environment. 

The Pierre Auger Observatory [SH] should be able to study the spectrum and 
composition of the ankle, a possible GZK feature in the UHECR spectrum and identify 
point sources of UHECR. The combination of any cosmogenic neutrino detection 
with detailed UHECR spectrum, composition and anisotropy analysis from the Auger 
Observatory would greatly help solve the mystery of the origin of UHECRs. 

6. Conclusion 

We studied the cosmogenic neutrinos production using a mixed composition hypothesis 
and several source evolution hypotheses. At high energies, the cosmogenic neutrino 
flux for our mixed composition model is very close to the pure proton case. In this 
range, ANITA and Auger can detect or strongly constrain the source distribution 
evolution with redshift |42| independent of the composition model assumed. In the 
lower energy range, the flux are higher due to interaction with the IR/Opt/UV 
backgrounds [T2|. The effect is strongest for pure protons, due to steeper spectral 
index, but it is also present in the mixed composition case mainly due to the direct 
photopion production process. If the confinement at UHECR sources is efficient, 
stronger fluxes can be expected in the Ice Cube energy range. In sum, Ice Cube, Auger, 
and ANITA will either detect cosmogenic neutrinos or provide important constraints 
on the confinement and the environment of the sources of extragalactic cosmic rays 
once combined with UHECR composition measurements. The combination of UHE 
neutrino and cosmic ray detectors will open a new window into the highest energy 
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phenomena of the present universe. 
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